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A canonical structure compatible with the action of the Lorentz group can 
be obtained considering the energy and time as conjugate variables of an 
extended phase space. Scalar probability waves, describing free relativis- 
tic particles, are associated with functional coherent states for an extended 
Liouville equation. Relativistic action waves are provided by distributions 
localized in the momentum space, evolving according to the continuity and 
Hamilton- Jacobi equations. Presuming the existence of minimum space and 
time intervals, the action distributions take the form of relativistic Wigner 
functions. The nonrelativistic quantum dynamics is retrieved approximating 
the time distribution function by a Gaussian wave packet. 
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1 Introduction 



The first evidence on the granular structure of the phase-space [1] emerged 
from the study of the relativistic system represented by thermal radiation, 
within classical statistical mechanics. However, the usual Hamiltonian dy- 
namics, of time-dependent coordinates and momenta, is not appropriate as 
a starting point for a Lorentz-covariant theory. 

A canonical structure compatible with the action of the Lorentz group can 
be obtained extending the nonrelativistic phase space by energy and time, 
as conjugate variables [2], setting thereby the framework to relate electro- 
dynamics, relativity, and quantum mechanics [3]. These new variables have 
been considered also in a variational approach to the coupled Vlasov- Maxwell 
equations [4]. In this work, relativistic probability waves are associated to 
the coherent solutions of the Liouville equation in the extended phase space. 

The extended Hamiltonian dynamics of a classical particle is reviewed 
in Sect. 2. As the observable time becomes a canonical coordinate, the 
evolution is described in terms of a parameter called universal time. The 
relativistic Liouville equation for the distribution function is presented in 
Sect. 3. Similarly to the nonrelativistic treatment [5], "action waves" are 
related to coherent functionals localized in the momentum space. For such 
probability distributions the Liouville equation reduces to the coupled con- 
tinuity and Hamilton-Jacobi equations. It is shown that for a distribution 
localized in a finite domain of the extended phase-space, the universal time 
is the expectation value of the time coordinate in the "intrinsic frame". 

The transition from action distributions in the extended phase space, to 
Wigner functions, arising from a discontinuous character of the inertial mo- 
tion, is discussed in Sect. 4. Following [5], discretization is introduced by 
presuming the existence of a minimum length £. Though, unlike the funda- 
mental length ip at the Planck scale (~ 10~^^ m) used in string theory [6], £ 
depends on the inertial parameter. The quantum "wave function" defined by 
discretization belongs to the extended Hilbert space presented in [7] , and with 
respect to the universal time, it evolves according to a relativistic Schrodinger 
equation. In the "stationary" case this reduces to the Klein-Gordon equa- 
tion, while the nonrelativistic limit is essentially " nonstationary" . Within 
the present approach, the Wigner transform for relativistic quantum systems 
can be defined directly as a quasiprobability over the extended phase-space, 
rather than over trajectories [8]. Conclusions are summarized in Sect. 5. 
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2 Classical dynamics in the extended phase 
space 

The phase-space M of a classical system can be extended to a phase-space 
M^, which includes the energy and time as conjugate variables [2]. The 
canonical coordinates on consist of the canonical coordinates on M, de- 
noted {qi,pi,i = l,n}, and {qo,po), supposed to be hnear functions of time 
and energy, qo = ct, respectively po = —E/c, where c is a dimensional con- 
stant, identified with the speed of light in vacuum [3]. 

Let u be the "universal time" parameter along the trajectories on M^, 
du = d/du the derivative with respect to u, and 

n 

Ch" = ^{duqi)dqi {duPi)dpi + {dj)dt {duE)dE (1) 

i=l 

the Lie derivative Cu-^f = fY- Here {*, *Y ^''^^ the extended 

Poisson bracket and Hamilton function, respectively, on M^. In the case of 
a nonrelativistic system can be taken of the form 

H%^H + cpo , (2) 

where H is the usual Hamilton function defined on M. For this expression, 
the corresponding equations of motion in the extended phase-space are 

dH , dH 

duQi = , duPi = (3) 
dpi dqi 

dH 

dj^l , duE^— . (4) 

The first group of equations reproduces the usual Hamilton equations on M. 
The second group shows that the choice of corresponds to u = t, and 
ensures the conservation of the energy when H is independent of time. 

In the extended phase-space, the transition from Newtonian to relativistic 
mechanics (recalled in Appendix 1) consists essentially in a change of Hamil- 
tonian. A free relativistic particle^ having q^ = {qo, q) and p^ = {po, p) as 
phase-space coordinates on = R^, can be described by Hq — —c^Pq — p^. 



^To include an external potential V(q) we may consider H'' = —c\J {po + V/c)^ — p' 
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When Pq ^ mgC^ ^ this extended Hamihonian reduces to the nonrela- 
tivistic expression i/^Q = p^/2mo + cpo, while in general it provides the 
equations of motion 

duQa = -c , = , duPo = (5) 

'Po - 

duq = c ^ , 4p = . (6) 
VPo - P^ 

With respect to a particular inertial frame, the usual velocity^ v = dq/ dt is 
the ratio v = cduq/d^qo = — cp/po- 

The inertial parameters for H^, defined by 

r^-^ , = 0,1,2,3 (7) 
If^ P^ op^ 

take the values 

I,^I^^I,^-Io^-^^mo , (8) 
as provided by the invariant value — moc^ of H^. 



3 Relativist ic action waves 

The statistical properties of classical systems composed of N identical rel- 
ativistic particles can be described by a distribution function > 0, de- 
pending on u, defined on the one-particle extended phase-space M*^. If ci?fi„je 
denotes the volume element around the point m'^ G Af^, then f'^{rrf',u) is 
normalized using the integrality condition 

j dVLm^r{rrf,u) , N>1 . (9) 

For macroscopic systems, the probability to find a particle localized in dQm'^ , 
proportional to f'^{m'^ ,u)6Qrn'' , is given essentially by the particle density in 
dfljn", while at small N we can expect a definition in terms of the average 
universal time interval of localization in dJl^e . 

^The usual Lagrangian on of the relativistic free particle can be found in [9]. 
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Let us consider a system consisting of a single relativistic particle (A^ = 1), 
with dQrn<= = d'^qd^p, {d^p = dpod^p, d^q = dq^d^q), and Hamiltonian 



Hq{p^) = —c\Jpq — v"^. The distribution function f^{q^.,p^.,u) evolves ac- 
cording to the relativistic Liouville equation 

dur+CH^,r = Q (10) 

where Ch^ provided by (1-6) has the form 

cp-V cpodo 

\IpI - VPo - 

with do = d/dqQ and V = d/dc{. Thus, (10) becomes 



\IpI - P^d^r - cpodor + cp-vr^o . (12) 

To find coherent solutions of this equation it is convenient to use the Fourier 
transform 

r (g^ k^, u) = j d'p e^^°«+^''■P/^(g^p^ u) , (13) 

which provides the spectrum of the "momentum frequencies", /c^/27r, = 
{ko, k). If ~ P^ can be expressed in the form moc + 5moc, where 6mo 
as a function of pi and is a power series, then by Fourier transform (12) 
becomes 

H'odur - ic^dkM' + ^cVk • Vf = , (14) 
where d^Q = d/dko, Vk = d/d]s., and formally 

H^, = -c^/i^tdij^'^i^tV^ . (15) 

Various densities in space-time, such as the localization probability n'^, or 
current J^, can be expressed directly in terms of / and its partial derivatives 
dkt, = d/dkn, /X = 0, 1, 2, 3, at /c^ = by 

n%q^,u) = ld'pnq^,p^,u)^nq^Au) , (16) 
J,{q',u) = — f d^pp,r{q^,p^,u) - -^d,,r{q',0,u) . (17) 



In general, the mean value of an observable 0{q'^,p'^), which is a polynomial 
as a function of the momentum components, has the expression 

{0){u) = J d^d^p Or{u) = J A -idke)r{q', O, u) . (18) 

A particular class of coherent solutions for the relativistic Liouville equa- 
tion (12) consists of the "action distributions" 

/o1?^P^^) = n%q^,u)5(po - doS)5(p - VS) , (19) 

where n'' is the localization probability density in space-time, and S{q'^,u) is 
the generating function of the Hamilton- Jacobi theory. By Fourier transform 
(19) becomes 

= ^egifcoaoS+ik.V5 ^20) 

while (14) reduces to the system of equations 



d^[n'^J{doS)^ - (V-S)2 ] = cdo{n'doS) - cV • (n'VS) (21) 

and 

n'd^J^O , /x = 0,l,2,3 (22) 

where = d/dq^, J ^ d^S - c^(doS)^ - (VS)^. Thus, the solution J = 
is nothig but the Hamilton-Jacobi equation in the extended phase-space, 

duS + H^{doS,VS) = . (23) 
Considering duS — moc^, (23) takes the form 

{doSy-{VSy = my , (24) 
and (21) becomes the continuity equation 

moduu'' = do{n''doS) - V ■ (n^VS) . (25) 
In terms of the density (16), the mean value of the time coordinate is 

{t) = ll d'q qo , (26) 
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so that 

du{t) = \j d\ qo dun' . (27) 

Let us presume that is hmitcd in time, confined to a finite volume V in 
space, and vanishes along the normal to the boundary of V. In this case 
(25,27) yield 

du{t) = — I dqJ d'q n'doS = (28) 

where {E) is the mean value of the energy. Because {E) is a positive 
constant^, (28) shows that {t) and u are in the linear relationship 

{t) = . (29) 

The localization to a finite domain in space-time makes possible to define 
(up to a translation) an "intrinsic frame" (IF), as the frame selected by 
the condition (p)iF = 0. Expressed in terms of the intrinsic expectation 
values, (29) shows that the universal time u corresponds up to the factor 
rriQC^ / {E)iY to the mean time in the intrinsic frame, (t)iF- Thus, if we can 
define tjiq as {E)ip/c^, then u = {t)ip. For a density n'^{q'^,u) = 6{qo — 
cu)n{q,u), localized in time, (25) reduces in the nonrelativistic limit to the 
usual continuity equation 5{t — u)[modun + V • (nVS)] — 0. 



4 The relativistic Schrodinger equation 

When the partial derivatives 9^5*, fi = 0,1,2,3, in (20) are written as fi- 
nite differences [S{q^ + 1^/2, u) — S{q^ — ^^/2, u)]/£^, the action distribution 
becomes 

/o^(g^r,«)= lim^(g'^ + ^,z.)**(g'^-^,«) , (30) 

where cr^ denotes the ratio = ifj./k^, and ^(g'^, it) is the complex function 
= \^exjp{iS/af^). However, if the limit of when both ifj, and /c^ 
decrease to zero is finite, having the same value a for all components, then 
we may consider also "quantum distributions" 

fl,(q',k^,u)^^(q'^ + ^,umq^-^,u) , (31) 

^As required to ensure the Lorentz action (61) [3]. 
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as possible functional coherent states for (14). In this case, the normalization 
condition (9) for the corresponding phase-space distribution /| takes the 
form 

I d'qd'pfUq^P^u) = I = 1 (32) 

and the phase-space overlap between two distributions /|,^, /^^ is 

< fUk / d'qd'pfUk = ^l^l^ (33) 

where 

(*l(^^l)|*2(^i2)) = /o^V^(g^^il)*2(g^^^2) . (34) 

A linear relationship = crA;^ with a finite, isotropic, Lorentz- invariant, 
universal phase-space clement a could be related in principle to the existence 
of minimum space and time intervals for a certain energy domain, but is more 
difficult to justify than in the nonrelativistic case [5]. In electrodynamics we 
can find limits such as the classical electron radius rg = ah/nieC = 2.8 fm, 
and the related cutoff energy e^/Aireore — nieC^- In general, we can note that 
in the intrinsic frame, after a cutoff at ~ Smgc^ of the energy range (Appendix 
2), still remains unchanged in 94% of the velocity domain [0,c). With 
this approximation, the inverse of (13) can be replaced by a multiple Fourier 
series in which the components of from the factor exp{—ikoPQ — i\<.-p) take 
an infinite set of discrete values separated hy k — 7r/3moC. Also, if the time 
distribution has the variance St"^, then the intervals of ordered, physical time 
(e.g. the lifetime tl = h/T for unstable particles) are greater than Sto, and 
the length between any two fixed endpoints along a trajectory parametrized 
by (t), greater than £ = c6to- Thus, a finite value a = £/k, nioc^Sto should 
be expected. For a quantum particle, with cr = ^ we get 5to ~ h/rriQC^ and 
£ ~ h/rrioc, proportional to the inverse of the mass. 

It is interesting to remark that beside the formal arguments, evidence for 
the physical relevance of the interval h/rriQC^ arises from the particle data. 
The values obtained for the ratio moC^/T ~ TL/Sto between the mass (in 
MeV) and decay width (F), using the experimental data [10] for meson and 
baryon resonances are represented in Figure 1, (A) and (B), respectively. 
These values are well interpolated by functions of the form 2.1 + C/F, where 
C is 1222 MeV for mesons and 1487 MeV for baryons. Thus, tl appears to 
be limited below by 2h/mQC^. 
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Denoting /| = {U^){U-^'i>*), with U = exp[a{kQdQ + k • V)/2], (14) 
becomes 

H^dJI, = i—[{um){U-'^*) - {U^){U-'m*)] , (35) 

where □ = — V^. 

A "static" distribution duf^ = is obtained if = a^, where a is a real 
constant. This constant can be estimated by using the expectation values of 
Hq or {HqY. For simphcity, {{H^Y) — rn^c^ means 

j d'qd'p {pI - - ml^)fl{q\ p^ u) = (36) 

/ A[(^*)** + *(/C**)] = , (37) 

where /C = — a^Q — ttIqC^. When U^i = a^, (37) yields aa^ = — mgC^, so 
that iC^ = 0, or 

-cr^n^^mgc^* . (38) 

In the quantum theory this represents the Klein-Gordon equation [11]. Al- 
though all particles described by (38) are unstable, closer to stability are 
quark-antiquark systems like the vr^ and mesons^ with a lifetime ~ 10 
ns, much larger than h/moc^ ~ 10^^^ s. 
In the nonstationary case (35) becomes 

H^[{Uid^^){U-^'^*) - {m)U-\idu^y)] = (39) 
[([/□*) (C/-^**) - (C/*)(C/-^n**)] , 

or 

//o'**_i9„*+ = , (40) 

//o'*+i9„*l = , (41) 

where ^+ = U"^ and ^'l = By contrast to the nonrelativistic case 

[5], in general (40-41) cannot be reduced to separate equations for ^+ and 

"^The singlet (triplet) states of e~e~^ positronium have a lifetime of 1.2 ns (140 ns). 
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^1 due to Hq, which acts on both functions. However, ^+ and become 
complex conjugate at A;"^ = 0, so that when a = h, the hmit 

lira {^*_)-^ H^^*_Uidu'^ ^ , (42) 

can be formally considered as a relativistic Schrodinger equation for the wave 
function Nonstationary solutions of this equation correspond for instance 
to wave-packets of the form 

'^{q'',u) = x{qo,u)ij{q,u) (43) 

where x{Qot'^) = XQo,PoiQo,u) is a Glauber coherent state [7] 

V Ca/tT 

with the centroid at Qq = —uVo/mo and variance c^/2f2^. The parameters 
Vo, Qo are related to the energy and time expectation values (E), (t) by the 
relations Vq = (po) = —{E)/c, respectively Qo = (qo) — c{t). 

The function f^{q^,p^,u) defined by f^{q^,k^,u) inveting (13) is^ 

rM.P^u) = — e-^'(^°-2°)'/^'-<='(^'°-^°)'/"''^V^(q,P,«) (45) 

where /^(q, p, ti) is the usual Wigner transform of ip{c\,u). Thus, the time 
variance 5t^ = (t^) — {ty = 1/20^, as well as the energy variance SE^ — c^Spl, 
Spl = {pD — {poY = o-^n^/2c^, are both finite, and satisfy the uncertainty 
relation 5E5t = a/2. 

With (44), the dependence on /cq in f^i^cf, k'^, u) can be separated in 

^o(A;o) = e-^^o*^o/2+iPofeo ^ (46) 

so that 

/|(g^^^^^) = ^o(A;o)lx(?o,^^)^(C>kV')(C>kV) , (47) 
and (40) becomes 

i/o'.5^feX*C/ki9„* = -— J^feX^C^kD* . (48) 
^Because /f remains finite for \pq\ < moc, (44) should be regarded as an approximation. 
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Here f/k = e'^'''^/^ and JF^ denotes the function J?-/,.(q, u) = go{ko)U^^iJj*{q, u). 

For nonrelativistic particles we can expect that ip evolves over a time-scale 
much larger than 1/ft, and approximate solutions of (48) can be obtained by 
taking the average over ^o- Using the equalities — (dux)'^ + xi^ui^), 

j dqoX*{qo,u)duX{qo,u) = "^'^"Qo = ' 

and 

(Po) = J dqoX*iqo, u)i-a'd'o)xiqo, u) = + 5pl , (49) 
the integration over qq in both sides of (48) yields 

H^o^kMi^du^ - = "^^FkMivD + a^V^)^' ■ (50) 

In general, H^J-'kUk is a complicated operator because Hq of (15) introduces 
mixed partial derivatives, acting both on ^o(^o) and C/k. However, in the 
limit A;*^ — > 



i^o ^o(A:o) ~ -c^(p2) + V2 ^o(A;o) 

where, according to the constraint (36), (pi) = mlc^ + (p^). Moreover, 
because 

hm I d'q{{p') + VDTkUx.^ - (p2) + £q ^V'^ = , 
we approximate H^J-'kUk ~ —moc^J-'kUk, so that when A;^ — >■ (50) becomes 

Using (49), this reduces further to 

ia^^^^--^{6pl + a'V')^|; . (52) 

^777.0 

The term Spl/2mo = ((jQ)^/4moc^ can be included in a global li-dependent 
phase-factor of ip, while by changing the parametrization to the mean time 
(t) — Qq/c, one obtains 

i(^d<^t)i^ = 7^^'^''^ ■ (53) 
^/ 
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Here, according to (49), Vq = —ymlc^ + (p^) with rrix = y^o ~ ^Po/c^ the 
effective mass, so that 

In the case of an atomic electron {a = h, rrix = rrie, 1 a.u. = a^mgC^), 
the correction term (He) = /i^(p^) (V^)/4mfc^ can be compared with the 
usual contribution due to the variation of mass with velocity, {H[) = —ip'^)/ 
8mgC^ [12]. For the ground state of hydrogen, (He) — —a'^/4: a.u., while 
{H[) — — 5q;^/8 a.u.. The whole correction to this order found by expanding 
in powers of a the exact solution of the Dirac equation is — a^/S a.u. [12]. 

The interval 2St = ^/(] = h/SE, is a measure of the time shift \Qq — 
Qo\/c for which the overlap \{xq'q,Vo\XQo,Vo)\'^ between two states (44), and 
the corresponding transition amplitude (34), remain significant. In general, 
this is much larger than 5to ~ h/rriQC^. For instance, if we take 5E fa eririxC^, 
where = Smx/rrix = 0.3 ■ lO^*' is the relative standard deviation at the 
measurement of the electron and proton mass, then 5t ~ 1.6 ■ 10^ 6tQ. In 
the case of the electron, Sto = i/c ^ 10~^^ s is comparable to the estimates 
of the "jump time" tj ~ 10^^° s for atomic transitions [13]. These change 
however the electron wave function over a distance larger than 10^£, so that 
5t could be a reasonable upper limit for tj. 



5 Summary and conclusions 

The phase-space description of the physical states provides the conceptual 
framework for nonrelativistic many-body theory, statistical mechanics and 
canonical quantization. The asymmetry between time and the usual phase- 
space coordinates requires though "the second quantization", to obtain a 
Lorentz-covariant quantum theory. 

For classical relativistic systems, we may also extended the usual phase- 
space by energy and time, as canonical variables. In this work, scalar proba- 
bility waves, describing free particles, are associated with functional coherent 
states for the Liouville equation (12) in the extended phase-space. 

The canonical equations of motion for a relativistic particle are presented 
in Sect. 2. As the usual time becomes a coordinate, the trajectories are pa- 
rameterized by a variable u called universal time. Action waves n^{q^,uY^\ 
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associated with specific coherent solutions (19) of (12), are discussed in Sect. 
3. These solutions are localized in the momentum space, and propagate 
according to the continuity (21) and Hamilton- Jacobi (23) relativistic equa- 
tions. It is shown that in a finite system the universal time is the expectation 
value of the time coordinate in the intrinsic frame. 

The transition from n'^(g^, u)^^^ to the quantum waves ^(g*^, u) is discussed 
in Sect. 4. Presuming the existence of minimum space and time intervals, 
the action distributions (20) take the form of the relativistic quantum dis- 
tributions (31). In the quantum case, the coherent solution of the Liouville 
equation is defined by the extended Wigner transform of the wave function 
provided by the relativistic Schrodinger equation (42). For an ideal non- 
interacting, "static" system, this reduces to the Klein-Gordon equation (38). 
Most physical situations are though nonstationary, as all mesons undergo 
irreversible decay, while in the nonrelativistic quantum theory time is the 
same as in classical mechanics. When time is described quasiclassically, by a 
coherent wave-packet, then over large intervals compared to the width, the 
extended formalism reduces to the usual nonrelativistic quantum dynamics. 

6 Appendix 1: Galilei and Lorentz actions 

Let us consider a particle with mass m, described by the Cartesian phase- 
space coordinates (q, p) . An infinitesimal Galilei transformation Fg : x 
R ^ X R, acting both on the coordinate space {Q — R^) and time (R), is 
defined by [q',t'] = [q,t] + 7(^, d, v, r)[q, t], where 



The algebra g of the Galilei group is isomorphic to so{3) + R^, and 7 e gf is 
specified by ^ e -50(3), d e R^, v e R^ and r e R. The parameters ^, d and 
V correspond to static rotations, translations and boost, respectively, of the 

space coordinates, while r describes translations along the time axis. 

The action Tq of the Galilei group can be lifted to an action F^ on the 
phase-space M = T*R^, by assuming that at the transformation specified by 
(55), the momentum also changes to 



7(C, d, V, r) [q, = [^q - d - tv, -r] 



(55) 



p' = p -I- — mv . 



(56) 
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However, as the boost transformations depend on time explicitly, and 



Po = Po + V • p/c 



(57) 



{po = —E/c), Tq can be lifted directly to an action Fm^ on the extended 
phase-space of Sect. 2. If the coordinates on are represented as 
column vectors 

Q 

then the infinitesimal transformation Tm<^ defined by (55), (56) and (57) takes 
the form of a canonical transformation 



q 




p 




, p = 









with 



p' 



X = 



p 



mv 




+ 



-Y 
-X 



+ 



C 

d 



Y = 



d 

T 



v/c 



p 







(58) 



(59) 



and 6 = c, 4 X 4 zero matrices. 

The mass m, introduced with the lift (56) is the positive, isotropic, inertial 
parameter for a Hamiltonian H defined on M. In the case of a Hamiltonian 

defined on M^, there is also an inertial parameter specified by the depen- 
dence of H'^ on the additional momentum component (po)- Following [3], this 
new inertial parameter is taken for simplicity as ±m, with -|- or — sign in the 
isotropic, respectively quasi-isotropic case. This yields a relationship of the 
form Po = ±mc, or E = ^mc^, which shows that the lift (58) of Tq should 
be obtained by placing the velocity v from (56) in the matrix a, instead of 
X. In this case, (59) is replaced by 



X 








Y 



d 

T 



v/c 



(60) 



According to (58), the new element in d changes the Galilei action (55) of 
the inertial equivalence group, and in the case > it provides the Lorentz 
action 

7z.(e,d,v,r)[q,t] = [eq-d-tv,-vq/c2-r] . (61) 

To find the action of the Lorentz group, (61) should be integrated to finite 
transformations. Let us presume that { = 0, d = 0, r = 0, and decompose 
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the vectors q, p with respect to the versor n of the boost velocity as q = 
q_L + gyn, p = p_L +P||n. In the representation 



we get a — pao, where p 







P± 


i\\ 




P\\ 






_ Po 


v|/c, 














(62) 



Ox = are 2x2 zero matrices, and 




1 



1 




(63) 



Because 



gPo-x _ ^.Qgjj p i + sinh p ax , 
then for a boost transformation with finite velocity V = Vn, the equations 



dq ^j,^ dp 
dp '^^ ^ dp 



aop 



(64) 



can be integrated to q'^ = q±, p'j^ = p±, and 

gy = cosh p gy — sinh p Qq , = cosh p qo — sinh p q\\ 

pjl = cosh p p|| + sinh p pq , Pq — cosh p po + sinh p py . 

These expressions show clearly the invariance of the Poisson bracket in the 
extended phase-space, because if {q^,quY — {Pf^^Pi^Y — {Qfi^PvY — ^ij.v^ 
then also {g^, q',Y = {p'^.p'^Y = 0, {q'^.p'^Y = 5^., Iim = 0, 1, 2, 3. 

The parameter p is related to the finite boost velocity V by physical 
considerations, such as y = cdq\\/dqo when dq'^^ — 0. The result V — ctanhp 
provides the standard Lorentz transformations 



qii - 



t' 



t-y- q/c^ 

1 - yVc2 



(65) 
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, Pll - VE E-V-p 
P I = ^^^^= > E' = . (66) 



For states with negative energy {E < 0), =Fv in (60) takes the — sign, the 
hyperbohc functions become trigonometric functions, and the Lorentz group 
SO (3,1) is replaced by the rotation group in space-time SO (4) [3], isomorphic 
to SU(2)xSU(2). 



7 Appendix 2: The relativist ic perfect gas 

For a nondegenerate gas of fermions with energy = ^J'^^~+rr^^ at equi- 
hbrium, the usual distribution function has the form [14] 

/m,t(p) = ^^e'-^-^^y'-^ , (67) 

so that if V denotes the confinement volume, then 

N^vJd'pf,,T{p) , E^vJd'pepf,,T{p) (68) 

arc the number of particles, and the total energy, respectively. The function 
(67) is also a stationary solution of the classical Fokker-Planck equation 

9t/ + -p-V/ = 7Vp-(- + A;BTVp)/ , (69) 
m m 

with m — Cp/c^. The energy can be expressed in the form 

where (?T(e) = e^^e^ — moC'^exp(— e/Zc^T). Here the upper integration limit 
Cm was presumed infinite, although in most physical situations particles with 
high enough energy can escape the system before thermalization. Moreover, 
the opening of pair creation reaction channels [15] at e = 3moC^, 5moC^, ... 
also affects the distribution. Therefore, a reasonable limit of the energy 
range for a stationary distribution with a well-defined number of particles is 
Em ~ SuiqC^, which corresponds to the maximum of gri^) at T — moC^/kB, 
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when the old sound velocity formula Vg = ykBT/rriQ yields Vg = c. 

By the simple limitation of the energy range, the inverse of (13) can be 
expressed in terms of a multiple Fourier series. Thus, a real function f{X) 
defined on the finite domanin [—A, A] can be represented in the form 

-1 oo 

m = ^ E , (71) 

n=— oo 

where 

fn = dX e'^^'^'^'fiX) . (72) 
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Figure 1. moc^/F from experimental data (*) and the interpolation functions 
2.1 + C/F (solid) for 32 meson resonances {to, rj, (f), tt, p, a, b, f) with F > 8.43 
MeV (A) and 45 baryon resonances (A/", A, A,S) with F > 15.6 MeV (B). 
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